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Abstract. System parameters should be tracked on-line to build a reconfigurable control systems even 
though there exists an abrupt change. For this purpose, a new performance index that we are studying 
is the speed of adaptation - how quickly does the system determine that a change has occurred? In 
this paper, a new, robust algorithm that is optimized to minimize the time delay in detecting a change 
for fixed false alarm probability is proposed. Simulation results for the aircraft lateral motion with 
a known or unknown change in control gain matrices, in the presence of doublet input, indicate that 
the algorithm works fairly well. One of its distinguishing properties is that detection delay of this 
algorithm is superior to that of Whiteness Test. 
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1. INTRODUCTION 

Designing a reconfigurable control system for an 
aerospace vehicle which operates throughout a 
wide range of flight environmental conditions that 
affect its dynamic characteristics is very attrac- 
tive. The problem of identifying aircraft system 
parameters from flight test data and their applica- 
tion to flight control has been successfully solved 
for linear models, by Balakrishnan (1972), Iliff 
(1973), and Maine (1986), when there are no se- 
vere changes in aircraft and environment. Even 
when the control design incorporates a degree of 
robustness, system parameters may drift enough 
to degrade its performance below an acceptable 
level. In other words, we still need to develop a 
theory for a better control system when aircraft 
flies at a high angle of attack where the linear 
model is no longer valid and when the robust con- 
troller does not work due to large variations in 
system parameters or noise levels. 

Adaptive controls are considered to be a promis- 
ing approach to give us a possible solution. An 
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adaptive control problem arises whenever system 
parameters are unknown or are subject to un- 
known variation which is usually a small change. 
It is also of increasing practical importance, since 
adaptive techniques are being used more and more 
in industrial control systems. However, this field 
is not mature yet, and only a set of scattered simu- 
lational results exists with no single quantitative 
measure for judging performances. Taking into 
account abrupt (or drastic) changes, which in- 
clude any faults in the system, in statistical mod- 
els appears as a natural complement of most of the 
adaptive techniques which track only slow varia- 
tions of parameters. So, the detection of abrupt 
change is essential in the design of a reconfigurable 
control system for these types of adaptive prob- 
lems. Therefore, the problem of detection, esti- 
mation, and diagnosis of changes in dynamical 
properties of signals or systems are considered, 
with particular emphasis on statistical methods 
for detection, to provide a general framework for 
change detection in signals and systems. 

The development of the change detection problem 
was stimulated by Wald (1947) when ‘Sequential 
Analysis’ was published and a Sequential Prob- 
ability Ratio Test was introduced. A ‘sequen- 



tial process inspection scheme’ was proposed by 
Page (1954) to detect a change in the mean by 
testing a weighted sum of the last few observa- 
tions, i.e. a moving average. Named the ‘Page- 
Hinkley’ rule, Page pointed out that this rule is 
equivalent to performing a SPRT. The theoreti- 
cal properties of this rule have been investigated 
for a long time from on-line and off-line points 
of view. The most significant works in that di- 
rection are Shiryayev (1978) and Lorden (1971). 
They formulated an ‘Optimal Stopping’ problem 
and proposed simple rules which are optimal in 
an appropriate sense. And the general problem 
of optimality of random processes in discrete time 
was considered and developed by Chow, Robbins, 
and Siegmund (1971). The problem of change 
detection has received growing attention during 
the last 20 years, as can be seen from survey pa- 
pers and books (Basseville, 1980; Basseville and 
Benvensite, 1988; Basseville and Nikiforov, 1993; 
Wilsky, 1976). But, the development of change 
detection method is still a relatively new subject. 

For reconfigurable control system, a new perfor- 
mance index that we are studying is the speed of 
adaptation - how quickly does the system deter- 
mine that a change has occurred? A new, robust 
algorithm developed in this paper is called Min- 
imal Time - Change Detection Algorithm (MT- 
CDA) which minimizes the time delay in detect- 
ing a change for fixed false alarm probability. 

A block diagram shown in Fig. 1 explains the re- 
configurable control system for the aircraft control 
that is equipped with a CDA. Before the CDA de- 
tects a change, the regulator is designed based 
on the initial parameter values. But after the 
change is detected at time r, the CDA block gen- 
erates an Activation Signal to re-initiate system 
identification. The system identification process 
finds a new set of parameter values proper for the 
changed system. The regulator uses the new pa- 
rameter values to build a new control law. This 
process is continued until the new parameter set 
converges. 


Algorithm complexity is definitely an issue for 
detection of each possible change and time of 
change. In spite of this complexity, we can have 
a control system which can detect its change and 
adapt it as quickly as possible even though the 
large amount of change is occurred at random 
time. Therefore, this MT-CDA will present a 
proper starting point for a reconfigurable control 
system in the theoretical point of view. 

2. MINIMAL TIME CHANGE 
DETECTION ALGORITHM 

2.1. MT-CDA: Detection of a Known Change 

The following state space model is considered. 

x n+ i = Ax n + ( Bi + (B 2 - Bi)x{n - 0))u„ 

v n = Cx n + ( Di + (D 2 - Di)x(n - 0))u n + GN n 

where x( TO ) = 1 if m > 0 and 0 otherwise. Based 
on this model, the basic problem to be solved is 
to detect a known change in control gain matri- 
ces ( B and D) as quickly as possible in order to 
build a reconfigurable control system whose pa- 
rameters could be switched from a known system 
to another known system after change detection. 
This is appropriate for a system which operates 
in two alternate driving points, and restricted to 
the class of models with a change of parameters 
which is known. The probability of making a false 
detection also needs minimized since the signals 
are disturbed by a white Gaussian noise. 

We want to formulate an algorithm which finds 
the detection time(r) to minimize the average 
time delay 

E[{r-e)+] ( 1 ) 

subject to a given false alarm probability, i.e. 

Pr.[r < 0] < a (2) 

where a is given. ‘Likelihood Ratio (change vs. 
no change)’, which we shall define to be the ratio 


Turbulence Nolee 



Joint density of «i, • • • , v n (with change) 

Measured output Joint density of iq , • ■ • , v n (without change) 

is introduced and p(-), Poo( ) denote the likeli- 
hood probability density function with and with- 
out changes at 9 — k. Then we define the follow- 
ing functional 


Ln,fc — 2 log 


p (vi,V2, - ■ ■ ,Vn\0 = k) 

Poo(vi,U 2 , • • •, V n ) 


( 3 ) 


l n - maxi,,*. 

k<n 


Fig. 1. A Block Diagram of Control System for 
Aircrafts with CDA 


and for each n 


( 4 ) 













Now a decision rule for CD A is 

r = min{n : £ n > h} (5) 

where r is a change detection time and h is a 
threshold value computed from the probability of 
false alarm equation (Eq.(2)). 

This algorithm stops the first time when l n ex- 
ceeds a given threshold. This is based on Eq.(3)- 
(5), so it can be regarded as ‘Maximum Likeli- 
hood’ treatment of the unknown change time, i.e. 
stop at t = n when for some k the observations 
v k> • • • > v„ are significant. 

By the contributions of Shiryayev (1978), Lorden 
(1971), and Kim (1993), there exists a change de- 
tection time and this procedure could be said opti- 
mal (at least asymptotically optimal), in the sense 
that it minimizes the detection delay for a fixed 
false alarm probability. 

2.2. MT-CDA: Detection of an Unknown Change 

In this section, minimal time detection of an un- 
known change instead of a known change is con- 
sidered, and its algorithm is derived. Multiple de- 
tectors are required to detect an unknown change, 
where the each detector is based on the known 
amount of change, in contrast to single detector 
required for the detection algorithm of a known 
change. The algorithm for the detection of only 
change time, which is a little bit simpler prob- 
lem, is first derived and the detection algorithm 
of both time and amount of change is followed. 

Detection of the Time of Change . The detection 
of change time is a significant problem in some 
applications, for example, 

a. Fault detection: all the processes are stopped 
when a failure occurs and are re-started after the 
failure is fixed. 

b. When a system requires the exact parameter 
set after a change, it is preferable to re-initiate the 
system identification just after detecting a change. 

In adaptive system sense, the 2nd case is a popu- 
lar one and a promising problem to be considered. 
To detect an unknown change, the CDA Block in 
Fig. 1 contains multiple change detectors. Each 
of change detection algorithm is based on a fixed 
known change, so each of detectors is basically 
the same one as derived in Section 2.1. By ex- 
amining the Activation Signal for each n, the real 
Activation Signal at time n is generated by the 
earliest one. The idea of this algorithm is that 
if there exists an unknown change, the exact (or 
the closest) detector will make Activation Signal 


at the earliest time. So the time of an unknown 
change will be detected with this algorithm if the 
amount of an unknown change is in the range of 
the multiple detector hypotheses. 

To formulate this MT-CDA for an unknown 
change, we consider the following system model 

z„ + i = Ax n + (By + (B 2 - Bi)x{n - 6))u n 

v n = Cx n + (D\ + (D 2 — Di)\(n — 9))u n + GN n 

where B 2 , and D 2 are unknown. The only differ- 
ence, compared to the Section 2.1, is that control 
gain matrices are unknown which stands for the 
unknown change. Like Section 2.1, we define the 
following functional, which is related to ‘Likeli- 
hood Ratio’, for each detector 

r i 2 , >ir' = k,B 2 = B' 2 ,D 2 = P 2 ) 

p(vi,v 2t - ■ ■ ,v n \0 > n) 

where each superscript ‘i’ stands for i-th detector 
(i = 1, 2, • • • , N) and for each n and each detector, 

4=ma xL' nk . (6) 

k<n * 

Then the decision rule for the detection of an un- 
known change is 

t* = min{n : C n > h 1 } (7) 

r = min r* (8) 

i 

where r* is a change detection time of each detec- 
tor, r is a real change detection time, and h' is a 
threshold value for each detector which is found 
from the probability of false alarm equation. So 
we can detect this kind of unknown change at 
time t which is the earliest detection time among 
all the multiple detectors considered in this al- 
gorithm. To have a good detection algorithm, 
the number of change detectors should be suffi- 
ciently large and the dynamic range of detectors’ 
change amount should be sufficient to cover the 
estimated range of changes in system parameters. 
This means that there exists a trade-off relation- 
ship between ‘Performance’ and ‘Complexity’ of 
the algorithm because the computational burden 
increases dramatically as the number of change 
detectors increases. 

Detection of the Time and Amount of Change . 

In some applications, it is insufficient to detect 
only the change time if the change is unknown, it 
is also required to detect the amount of change. 
Fig. 2 is a modified block diagram of Fig. 1 for 
this purpose. CDA Block is a change detection 
algorithm which can estimate the amount of an 
unknown change and the change time. It gener- 
ates the Activation Signal(A5) and an Estimated 



Degree of Change (EDC) which is closely related 
to the real amount of change limited by the reso- 
lution of the detectors. When the system change 
is detected, the system parameter set is updated 
based on EDC and Parameter Update Block that 
rules the parameter update law. There is no 
necessity of System Identification routine. This 
looks simpler than the adaptive control system 
explained earlier, but the CDA Block is much 
more complicated than CDA Block. So, Fig. 2 
could be considered another version of block di- 
agram for the total adaptive control system that 
is appropriate for a system that undergoes an un- 
known change in parameters. To determine the 



Measured Output 


Fig. 2. A Block Diagram of Control System for 
Aircrafts with CDA 


amount of change, a ‘Decision Mechanism’ is re- 
quired. Based on the A5, each of that comes from 
a known change detector, the Decision Mechanism 
generates a real one and an EDC at time n. The 
algorithm formulation is the exactly same as be- 
fore (Eq.(6)-(8)) except the Decision Mechanism, 
so only the details of Decision Mechanism is ex- 
plained. 


min{n : t' n > h'} 

( 9 ) 

minr' 

i 

(10) 

min{i : r 1 } 

(11) 


From Eq.(9)-(11), an unknown change could be 
detected at time r whose EDC is equivalent to 
the K-th detector, if the amount of an unknown 
change is in the range of the amount of multiple 
detectors. The difficulty of the decision of the ex- 
act EDC exists in this decision rule, called D.R. 

# 1, because some of change detectors could re- 
spond at the same time. That means the necessity 
of a better decision rule . Since the exact value 
of i l n is rapidly increasing, the following new de- 
cision rule (D.R. # 2) which is the modification 
of Eq. (11) could be proposed. 

f max{i : (^ — h*)} if 3 more than 2 min r* 
K ~ ( min{i : r*} otherwise 

There exists the same trade-off between perfor- 
mance and the complexity of the algorithm. Only 


a finite number of change detectors are used in 
real simulations. 


3. SIMULATIONS AND RESULTS 

A lateral rigid body dynamics of PA-30, which is 
a Linear Time-Invariant, is used for simulations 
(Kim, 1993). 

3.1. Detection of a Known Change 


The mean(M[-]) and the variance(V[-]) of detec- 
tion delay (r — 9), with respect to the different 
change time(0), based on 100 runs hereafter, are 
summarized in Table 1. 50% change is detected 
earlier than 25% or 10% change. From the re- 
sult of Table 1, it is shown that different known 
changes could be detected with some delays. 


Table 1 Results for Various Change Times 
when Change is in Rolling Motion 


9 

10 

20 

30 

40 

50 

60 

M[t — 6] 

7.6 

0.9 

7.1 

4.0 

11.7 

10.3 

V[t-0] 

0.4 

0.3 

4.5 

0.4 

26.6 

8.2 

(a) 10% Change 

9 

10 

20 

30 

40 

50 

60 

M[t - 9} 

6.2 

0.1 

3.1 

3.0 

1.6 

4.6 

V[r - 9] 

0.5 

0.1 

0.6 

0.1 

1.9 

1.0 

(b) 25% Change 

9 

10 

20 

30 

40 

50 

60 

M[t - 9} 

4.5 

0.0 

1.4 

2.2 

0.2 

2.3 

V[t-9] 

0.7 

0.0 

0.3 

0.2 

0.1 

0.5 


(c) 50% Change 


3.2. Detection of an Unknown Change: Time 


Simulations with 6 change detectors are done first 
and then combined to analyze their results as 
summarized in Table 2. From this Table, the fol- 
lowings are observed: 

• For a given detector, large amount of change 
is detected earlier than small amount of 
change. 

• Rolling motion: 

- 10% Change: 50%D > 25%D > 10%D 

- 25% Change: 10%D,50%D > 25%D => 
10%D,50%D ~ 25%D 

- 50% Change: 10%D > 25%D > 50%D 
=> 10%D ~ 25%D ~ 50%D 










Table 2 Summarized Simulation Results for 9 = 30 


where R.D. 

, Y.D., R.C. 

Y.C., 

and N.D. stand for 

Rolling Detector, Yawing Detector, Rolling Change 

Yawing Change, and No Detection, respectively 



10% 

25% 

50% 


R.D. 

Y.D. 

R.D. 

Y.D. 

R.D. 

Y.D. 

io% R.C. 

7.05 

16.2 

9.02 

15.8 

12.3 

15.7 

io% Y.C. 

N.D. 

16.6 

N.D. 

16.6 

N.D. 

16.5 

20% R.C. 

3.88 

15.6 

3.84 

15.6 

4.31 

15.0 

20% Y.C. 

N.D. 

16.6 

N.D. 

16.6 

N.D. 

16.5 

25% R.C. 

3.23 

15.2 

3.11 

15.0 

3.29 

14.9 

25% Y.C. 

N.D. 

16.6 

N.D. 

16.6 

N.D. 

16.5 

30% R.C. 

2.73 

15.0 

2.55 

14.9 

2.73 

15.8 

30% Y.C. 

N.D. 

16.6 

N.D. 

16.6 

N.D. 

16.5 

40% R.C. 

2.18 

14.8 

1.90 

14.7 

1.94 

14.8 

40% Y.C. 

N.D. 

16.6 

N.D. 

16.6 

N.D. 

16.5 

50% R.C. 

1.81 

14.7 

1.57 

14.7 

1.43 

14.7 

50% Y.C. 

N.D. 

16.6 

N.D. 

16.6 

N.D. 

16.5 

60% R.C. 

1.55 

14.5 

1.25 

14.6 

1.11 

14.6 

60% Y.C. 

N.D. 

16.6 

N.D. 

16.5 

N.D. 

16.5 


where ‘D’ denotes a detector and the detec- 
tion delay is quite reduced when 10% detec- 
tor is used to detect 10% change but there 
is no big difference in detection delay of 
3 change detectors to detect 25% or 50% 
change. 

• Yawing motion: The results are the same as 
rolling motion, but detection delay is large 
and the difference of delay with respect to 
each change detector is small since the por- 
tion of yawing motion in PA-30 is small. 

From the above observations and the trade-off be- 
tween the number of change detectors and the 
computational burden, it could be concluded that 
10% change detector is sufficient when the detec- 
tion of change larger than 10% is concerned. 

3.3. Detection of an Unknown Change: Time 
and Amount 

Simulation results, with D.R. # 2, for 0 = 40 are 
summarized in Table 3. From this table, it could 
be shown that the unknown changes at 9 = 40 are 
fairly well detected except 30% and 40% changes, 
which might be fixed if the better decision rule is 
applied. Simulation results for 6 = 10, 20, and 30 
have the same characteristics. Therefore, it could 
be concluded that the given unknown changes are 
detected in statistical sense, though there exists 
a difficulty in deciding the exact EDC in some 
situations which might be improved by applying 
the better decision rule. 

3.4. Change Detection via Whiteness Test 

To examine the validity of MT-CDA, simulations 


Table 3 Simulation Results for 9 = 40 



Detection Number 

M[t - B] 

V[r -9} 

10% D 

69 

4.12 

0.31 

20% D 

29 

3.55 

0.39 

30% D 

0 



40% D 

1 

4.0 

0.0 

50% D 

1 

2.0 

0.0 

(a) 10% change 


Detection Number 

M[t - 9] 

V[r - B] 

10% D 

0 



20% D 

87 

3.15 

0.15 

30% D 

3 

3.33 

0.22 

40% D 

3 

3.0 

0.0 

50% D 

7 

2.29 

0.49 

(b) 20% change 


Detection Number 

M[t - 0 ] 

V[t - 9] 

10% D 

0 



20% D 

18 

2.89 

0.10 

30% D 

27 

3.07 

0.07 

40% D 

33 

2.97 

0.03 

50% D 

22 

2.23 

0.18 

(c) 30% change 


Detection Number 

m[t - e\ 

V[t - 6] 

10% D 

0 



20% D 

4 

2.5 

0.25 

30% D 

3 

3.0 

0.0 

40% D 

23 

3.0 

0.0 

50% D 

70 

2.36 

0.23 

(d) 40% change 


Detection Number 

M[t - B ] 

V[t - B] 

10% D 

0 



20% D 

1 

2.0 

0.0 

30% D 

0 



40% D 

2 

2.0 

0.0 

50% D 

97 

2.23 

0.18 


(e) 50% change 


of Whiteness Test are done and results are com- 
pared in this section. Only the detection problem 
of a known change is considered here. 

Table 4 shows the detection delay of each algo- 
rithm for various change time($). For every tests, 
the detection delay of MT-CDA is smaller than 
that of Whiteness Test. Therefore, the simulation 
results show that MT-CDA is superior to White- 
ness Test in the sense of detection delay compar- 
ison. 

4. CONCLUSIONS 

To design a reconfigurable control system, a sim- 
ple model with a single change in parameters of 
control gain matrices at random time whose a 
priori distribution is known has been considered. 
Changes in these matrices mean the failure in the 
actuators, so possible damages in the control sur- 
faces and/or one of the engine failure could be 



Table 4 Simulation Results for Various Change 
Times when 25% Change is in Rolling Motion 


0 

10 

20 

30 

40 

50 

M[t - 6} 

6.2 

0.1 

3.1 

3.0 

1.6 

V[r - 0] 

0.5 

0.1 

0.6 

0.1 

1.9 

(a) Minimal Time-Change 

Detection Algorithm 

e 

10 

20 

30 

40 

50 

M[t - 0} 

11.7 

4.0 

17.9 

7.9 

37.4 

V[t - 0} 

0.2 

0.0 

0.1 

0.1 

4.8 


(b) Whiteness Test 


analyzed with this model. Therefore, consider- 
ing changes in the control gain matrices is the 
proper starting point before we consider all kinds 
of changes. MT-CDA, based on Likelihood Ratio 
Test, for this simple model has been developed 
whether a change is known or unknown and it has 
been shown to be optimal for a fixed false alarm 
probability. 

Simulation results for lateral motions show that 
a known change can be detected with some de- 
lays which should be minimal for the given con- 
dition as theory indicates. From simulation re- 
sults for the detection of change time, it could 
be concluded that 10% change detector, which 
is the smallest one, is sufficient and it is robust 
in the sense that it could detect any change(> 
10%) with sub-minimal detection delay. And the 
unknown change could be detected, in statistical 
sense, though there exists a difficulty to decide 
the exact EDC in some situations. The better 
decision rule needs to be developed to solve this 
difficulty. The algorithm with multiple change de- 
tectors, if implemented in full, will yield the best 
performance for the widest class of change. But 
full implementation is another problem. Simu- 
lations have been done throughout the various 
situations and the results are meaningful (Kim, 
1993). Detection via Whiteness Test has been 
also formulated for algorithm comparison in the 
sense of detection delay. MT-CDA is superior to 
the Whiteness Test in this sense. 

In this paper, MT-CDA has been developed and 
demonstrated by simulations when a change is 
in control gain matrices. With this MT-CDA, a 
reconfigurable control system which is a highly 
adaptive system shown in Fig. 1 or Fig. 2 could 
be built and work throughout a wide range of con- 
ditions. 
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